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ABSTRACT 

It is well known that collisionless shocks are major sites of particle acceler- 
ation in the Universe, but the details of the acceleration process are still not 
well understood. The particle acceleration rate, which can shed light on the ac- 
celeration process, is rarely measured in astrophysical environments. Here we 
use observations of gamma-ray burst afterglows, which are weakly magnetized 
P_f relativistic collisionless shocks in ion-electron plasma, to constrain the rate of 

particle acceleration in such shocks. We find, based on X-ray and GeV after- 
glows, an acceleration rate that is most likely very fast, approaching the Bohm 
limit, when the shock Lorentz factor is in the range of V ~ 10 — 100. In that case 
X-ray observations may be consistent with no amplification of the magnetic field 
in the shock upstream region. We examine the X-ray afterglow of GRB 060729, 
which is observed for 642 days showing a sharp decay in the flux starting about 
400 days after the burst, when the shock Lorentz factor is ~ 5. We find that 
inability to accelerate X-ray emitting electrons at late time provides a natural 
explanation for the sharp decay, and that also in that case acceleration must be 
rather fast, and cannot be more than a 100 times slower than the Bohm limit. 
We conclude that particle acceleration is most likely fast in GRB afterglows, at 
least as long as the blast wave is ultra-relativistic. 

Subject headings: 



Introduction 



Astrophysical collisionless shocks are efficient particle accelerators. The signature of 
ultra-relativistic particles that are accelerated in these shocks is seen in a variety of astro- 
physical phenomena and over a wide range of environments. Nevertheless, despite of an 
extensive study, the acceleration processes are still largel y unknown. One of the leading 



candidates is the diffusive shock acceleration (DSA; e.g., iBelll Il978l ; iBlandford &: Ostriker 
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19781 ; iBlandford fe Eichlerl 119871 ). where charged particles are accelerated by crossing the 
shock back and forth. The acceleration time in DSA depends on the duration that it takes 
a particle to close a single cycle, i.e., to cross the shock back and forth one time. This time 
increases with the particle energy, and it sets the maximal Lorentz factor, ^ max that a parti- 
cle can achieve. Thus, measuring ^ ma x provides a direct information about the acceleration 
process and about the physical conditions in the acceleration site. 

The reflection of particles back and forth through the shock is believed to be done by 
scattering on fluctuating magnetic fields. This process is typically approximated as a dif- 
fusion in direction of the particles velocity, in which case the duration of the acceleration 
depends on the particle mean free path, A. It is reasonable to assume that typically the 
shortest possible mean free path is the Larmor radius r\ (Bohm limit), and thus to param- 
eterize the diffusion by A = rjri where rj > 1 is generally expected, although this is not a 
hard lower limit. The value of rj, which measures the diffusion efficiency and how fast is 
the acceleration, was constrained only in a small number of systems. Probably the best 
estim ate is obtained from the gamma-ray spectrum of the Crab nebula (jde Jager fe Harding 
19921 ). The spectrum shows two components, where the lower energy component, which 
is most likely dominated by synchrotron, shows a cut-off around 100 MeV. The fact that 
synchrotron emission reaches these energies, given the rapid synchrotron cooling, implies 
rj « 1. H igher energy synch rotron emission during flares suggest that maybe even rj < 1 is 
required (lAbdo et al. Il201ll ). Another system w here r? was claimed to be measured is the su- 
pernova remnant SNR RXJ1713. 72-3946, where lUchiyama et al. I (120071 ) find that observed 
X-ray variability on a year time scale indicates on rj ~ 1. Thus, two very different acceler- 
ation sites, one relativistic, possibly highly magnetized shock in pair plasma and the other 
Newtonian shock in ion-electron plasma suggest that particle acceleration, if dominated by 
DSA, is extremely fast. Calculating 77 from first principles is impossible at this point, since 
it depends on the unknown shock structure, and most importantly on the properties of the 
upstream and downstream magnetic fields. Calculations of 7 7 in relativistic shocks were done 
only by assuming the mag netic field structure. For example, iLemoine fe Pelletierl (120031 ) and 
Lemoine fe Revenul ( 120061 ) find that when a Kolmogorov magnetic turbulence spectrum is 
assumed in the upstream region, then rj w 10. 

Here we examine the constraints that can be obtained on rj from observations of cut- 
off, or the lack thereof, in GRB afterglow light curves and spectra. These afterglows are 
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studies of particle acceleration in GRB afterglows assumed 77 = 1 — 10 and used the lack 
of spectral cut-off in observed X-ray and GeV emission to constrain the magnetic fi eld up- 



stream and/or downstream of the shock (ILi fe Waxmanll2006l ; iPiran fe Nakarll2010l [PN10]; 
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Barniol Duran fc Kumarll2011al ; iLi &: Zhaoll201ll ). Here we take a different approach asking 



how well 7] can be constraint. Moreover, we find that the fast decay observed in the extraor- 
dinarily year long X-ray afterglow of GRB 060729 is naturally explained by the inability of 
the shock to accelerate X-ray emitting electrons, providing a measurement of rj in that case. 

In section £j2]we describe the various limits on the ^ m ax- The resulting limits on t], for 
various circum burst density profiles, are derived in £j3l The special case of the GRB with 
the longest duration X-ray afterglow, GRB 060729, is discussed in §U 



2. Limits on the maximal Lorentz factor 

The two main factors that limit the accelera tion of particles in a decelerating relativis 



tic blast wave are confinement and cooling (e.g., ILi fc Waxmanll2006l ; iPiran fc Nakarl 12010 



Barniol Duran fc Kumar! l2011af ). Below we shortly discuss these limits (see PN for detail). 
Observations indicate that the magnetic field in the downstream region is amplified in GRB 
external shocks well beyond the effect of compression. Thus, most of the emission take place 
in the downstream region while a particle spends most of its acceleration time in the up- 
stream region. Throughout the paper we assume that rj is similar in the shock downstream 
and upstream, but we highlight which of the observations constrain r\ in downstream region 
and which constrain rj in the upstream region. 

Confinement is limited by the ability of the accelerated particle that is moving in the 
upstream region to cross the shock back into the downstream region. Thus, confinement is 
limited by r\ in the upstream region. Its limit on the maximal Lorentz factor is set by the 
requirement that the particle complete a turn of 180°, as seen in the shock frame, while the 
shock propagates a distance f u R, where f u accounts for the shock deceleration (see PN10). 
Thus, 

1'conf « -^y„* (1) 

1 r]m e c z 

where " ' " denotes quantities in the shock rest frame, R is the shock radius, e and m e are the 
electron charge and mass, c is the speed of light and B u is the rest frame upstream magnetic 
field. 

The two processes that dominate cooling are synchrotron and inverse Compton (IC). In 
both cases the maximal Lorentz factor is found by equating the acceleration time (i.e., the 
time to complete a Fermi cycle in relativistic shocks) to the relevant cooling time. In case 
of synchrotron: 

, / 6ne \ 5 
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where <jt is the Thomson cross-section and Bd is the downstream magnetic field. Since obser- 
vations indicate that the downstream magnetic field is amplified by the shock, synchrotron 
cooling in the downstream region is more limiting than in the upstream region and it sets a 
limit on rj in the downstream region. 

Inverse Compton cooling is more efficient in the upstream region, since the radiation 
field is similar in both sides of shock, but a particle spends more time in the upstream region, 
where the magnetic field is lower. Thus, 7j C is limited by 77 in the upstream region: 

v ~ ( 3eK V = v ( K V /2 ™ 

lic ~ WtUU< v'kn)) ~ %vnch \Y{iw)B d ) ' 1 ' 

where B' u is the magnetic field in the upstream region as measured in the shock frame (related 
to the rest frame magnetic field by B' u TB u ). U' rad {< v' KN ) is the shock frame radiation 
energy density at frequencies smaller than: 

JknWio) = (4) 
h lic 

where h is the Planck constant. Equation [3] gives also the relation between ^' IC and j' synch 
using the ratio between IC and synchrotron cooling rate in the downstream region, Y^Jjq) = 
UU< v' KN )/(Bj/8n). 



3. Maximal observed frequency and limits on rj 



Below we derive the constraints that X-ray and GeV afterglows set on 77. We assume that 
the observed emission is synchrotron radiation generated by a quasi-spherical decelerating 
adiabatic blast wave. This is almost certainly the case in many X-ray afterglows, at least 
during the first day. The origin of the observed long lasting GeV emission, which is seen up 
to ~ 1000 s after some bursts, is still unclear, althoug h observations suggest that it i s also 



synchrotron emission f rom the decelerating blast wave (IKumar fc Barniol Duranll2009l . 12010 



Ghisellini et al. Il2010l ). Since afterglow observations suggest that the circum-burst density 



profile vary from one GRB to another, we consider here two typical external density profiles, 
one constant, as expected for the interstellar medium (ISM) and one oc R~ 2 as expected for 
a stellar wind. 



3.1. ISM 



Under the assumption of spherical expansion in a constant density n, the radius and 
Lorentz factor of an adiabatic blast wave with energy E at an observer time t are (e.g., 



-5 - 



Sari. Piran fc Narayanlll998l ): 



i?«6-10 17 cm (^) 1/4 (i±i)- 1/4 ^ 4 

(5) 

r ^i2 (i^f* (i±*y/* t -v* 

where z is the burst redshift and q x denotes the value of q/10 x in c.g.s. units. We assume 
that the downstream magnetic field is a constant fraction, 634-, of the internal energy behind 
the shock so Bd ~ (32TTEB,dnr 2 m p c 2 ) 1 ^ 2 , where m p is the proton mass. If the magnetic field 
in the upstream is not amplified by a precursor to the shock then it is expected to be 
constant and of order of 10/iG. If it is amplified then it may be significantly larger. Finally, 
f u = 1/3 in ISM (PN10). Using the synchrotron emission from the downstream region, 
hv = T r y' 2 eB ( i/(27rm e c), we obtain the maximal frequency that is dictated by the limits 
discussed above: 

hv * « 2 • 10 10 eV n~ 2 £ 3/ V /2 B 2 rT 1/4 r 1/4 

hv sync « 2.5 • 10 8 eV 17-H5 3/8 (6) 

/m/jc ~ 1.5 • 10 5 eV Ty-^^-se-^^Y- 1 ^ 3 / 8 

where here, and throughout the paper, we derive values for the typical z — 2. We also ignore 
dependencies on parameters that are raised to the power of 1/8 since these cannot affect the 
result by an order of magnitude, which is the accuracy of our calculation to begin with. The 
constraints that we derive for confinement and synchrotron are similar to those of PN10, 
which assumed p = 1 , and the constraint on the IC cooling is similar to the one derived in 



Li &: Waxmanl (120061 ) . which used a canonical value of rj = 10. 



Equation [6] implies that for canonical GRB parameters confinement does not play an 
important role when the afterglow blast wave propagates into ISM, with the possible ex- 
ception of very early time GeV emission^. The limit provided by synchrotron cooling is the 
most robust as it is independent of almost anything^], except for the time since the explosion, 



^or early time GeV photons, the confinement limit can be of the same order of magnitude as the 
synchrotron cooling limit. However, it is less robust than synchrotron limit due to poorly constrained 
parameters such as B u . Note, that in specific cases, confinement may provide the most stringent constraint. 
An example is the extreme case where the m agnetic field is assumed not to be amplified by the shock and 
e B ,d ~ 10" 8 |Kumar fc Barniol Duranlbood PN10.) 

2 The synchrotron limit on the rest frame frequency is independent of anything and is in general hv' 
50MeV l|de Jager fc Hardindll99l lLvutikovll2009l: iKirk fc Revillel lioiol ) . The observed frame limit depends 
only on the Lorentz boost and cosmological redshift, which vary by less than an order of magnitude for 
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which is typically well measured: 



V<3 



hp, 



obs 



1 Gov 



100 s 



-3/8 



(7) 



Several Ge v photons are seen ~ 100 s after the burst sta rts in a number of the Fermi-LAT 
GRBs fe.g. JAbdo et al. ll2009a] Jbl: I Ackermann et al. 1120111 ) while > 100 MeV photons are seen 
in large numbers up to ~ 1000 s after the burst in many Fermi-LAT GRBs. Therefore, if the 
Gev emis sion is emitted by synchrotron proces s in the external shock, as suggested by several 
authors ( Kumar fc Barniol Duranll2009l |2010| ; iGhisellini et al. 1120101 ) then the acceleration 
process in ultra- relativistic (r ~ 100), weakly magnetized, shock must be extremely fastjfl 
with T] < 1 in the shock downstream. 

The IC constraint depends on the value of the Y parameter. Ya e v, the Y parameter 
of Ge V emitting electrons, varies by many orders of magnitude across the relevan t phase 
space (TNakar. Ando fc Sarilhood : iLi fc Zhaoll201l[ iBarniol Duran fc Kumarl ROllbh . Asa 
result the synchrotron cooling limit is more stringent in part of the phase space. Therefore, 
considering the robustness and tightness of the synchrotron constraint, we do not attempt 
to cover here the possible IC limits on the GeV emission. However, when considering X-ray 
emission the synchrotron limit is very loose. Therefore, we consider the IC limit in that case, 
for which we need to evaluate Y x . This is not trivial due to Klein-Nishina [KN] effects that 
play different roles over various areas of the phase space. An upper limit on Y x can be easily 
obtained by assuming that the electrons are in the fast cooli ng regime and tha t KN effects 
are negligible. In that case Y x = \/e e /eB,d if £ e > £34 (e.g., ISari fc Esinll200ll ). where e e is 
the fraction of the internal energy behind the shock that goes into accelerated electrons. For 
typical parameters the electrons are cooling slowly at t > 10 4 s, implying that Y x is smaller 
for electron distribution with a power-law index p > 2: Y x = y s e / £ B,d(lc/lm) 2 were 7 m is 
the typical (also minimal) Lorentz factor of accelerated electrons and 7 C is the Lorentz factor 
of electrons that are cooling over dynamical time scale. In addition, over a large range of 
the parameter space KN suppression can be important, reducin g the value of Y x further. T o 
account for these effects we use equations 46, 59, 60 and 63 of iNakar. Ando fc Saril ( 120091 ). 
which take consideration of the KN effects and their feedback on the electron distribution, 



reasonable values of E, n and z. A frequency that is higher by some factor may be obtained in configurations 
where the accel eration takes pl ace in a relatively weak magnetic field and then the radiation takes place in 
a stronger field rtLvutikovll2009h . In that case the synchrotron limit on rj is larger by that factor. 



(l2010l) . 



then its Lorentz factor 



3 Note that if the blast wave is radiative, as suggested by IGhisellini et al. 
drops faster with time than in the adiabatic case. Since the synchrotron limit on rj depends only on the 
Lorentz factor, the limits that it provides in a case of a radiative blast wave are tighter than those that we 
derive here for an adiabatic blast wave. 
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1 /2 

to scan the phase space for the value of s^ d Y x (which appears in the IC limit of equation 
|6j). We scan the parameter phase space and find that if the fraction of downstream region 
internal energy that goes t o electrons is e P = .1, the electron distribution power-law index 
is in the range p = 2 — 2.8 ( Curran et al. 2010 ) and Es,d > 10~ 3 , then the value of e^ 2 d _ 2 Y X 
is typically in the range of 0.3 — 3 and its dependence on the other parameters, n, E and t, 
is rather weak (most of the dependence in this range is on p due to the fraction of energy 
that is in fast cooling electrons, while KN suppression is rather mild). When eb <C 10 -3 

1 /2 

and/or e e <C 0.1 Klein-Nishina effects significantly suppress Y x and Eg Y x -C 1. We therefore 
conclude that for the canonical values of Eb > 10~ 3 , E e = 0.1, n w 1 and p = 2 — 2.8 



10 kev , 
V < 15 -r B 



",-5^5 



3/8 



in the shock upstream region. If, however, Eb "C 10 3 or e e C 0.1 or n <C 1, then X-ray 
emission does not provide strong constraints on r\. 

Many afterglows show X-ray emission (0.2-10 keV) that is brig ht for days and in 



some cases weeks, w ithout showing a clear sign of spectral softening ( ILiang et al. 1 12008 



Racusin et al. II2009I ). Thus, since afterglow modeling typically implies Eb > 10 3 or £ e w 0.1 



and an ISM circum burst environment, observations of X-ray afterglows suggest that the ac- 
celeration mancinism in relativistic shocks (r ~ 10 — 50) is fast. As evident from equation 
[HI this limit depends on various parameters. Some are constrained rather well, e.g., e e , while 
others are less constrained, e.g., no- Most important is the dependence on B u . There is a 
viable possibility that the interaction of accelerated particles that run ahead o f the shock 



significantly amplifies the u pstream magnetic field (IBlandford fc Eichlerl 119871 ; iBelll 12004 



Milosavljevic fc Nakarll2006l ). If this is the case th en the limits provided by X-ray observa- 



tions are rather loose. In fact. lLi fc Waxmanl (120061 ) concluded, based on X-ray observations, 
that the upstream magnetic field must be amplified at least up to > Oln^mG. This con- 
clusion was based on the assumption that rj = 10. They also take as a canonical value 

1 /2 I 

£ Bd-2*x ~ v 10, which they calculate by considering only the part of the phase space where 
KN effects are negligible and by taking p = 2, for which, Y x is not suppressed by the slow 

1/2 

cooling of m ost of the elect r ons. O ur results show that a more careful estimate of s B d _ 2 Y x 
reduces the iLi fc Waxmanl (120061 ) limit by at least a factor of a few. In addition, if the 
acceleration is as fast as suggested by the recently detected GeV emission, and rj ~ 1 also in 
the upstream region then the limits on the upstream field drop to [iG level, implying that 
current X-ray observations do not provide strong evidence for magnetic field amplification 
in GRB afterglows. 
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3.2. Wind 

The mass density profile in a wind from massive stars is p = AR~ 2 . Under the assump- 
tion of spherical expansion the rad ius and Lorentz facto r of an adiabatic blast wave with 



energy E at an observer time t are ([Chevalier Sz Lil I200CH ) : 

, 1/2 fl+*\-V2+l/2 

(9) 



J R = 3-10 17 cm(ff) 1/2 (i^)" 1/2 4 



r = i2 (^y\^) l,A t - 1/A 



where A* = 5 . 10 i^ g / cm ■ Similarly to the ISM case we assume that the downstream magnetic 
is a constant fraction, £34,-2 of the internal energy behind the shock. Contrary to the ISM 
case, the upstream magnetic field is not constant. The magnetic field in the upstream region 
depends on the wind magnetization and flux freezing implies B u oc -R -1 , assuming that 
upstream field is not amplified by the shock precursor. The normalizat ion depends on the 



wind velocity and on the surface rotation velocity and magnetic field ( jGoldreich &: Julian 



19701 ). which are not tightly constrained. For typical parameter s of a Wolf-Rayet wind a 



field of ~ 10/xG is expected at R = 10 18 cm (jEichler fc Usovlll993r ). but it can be more than 



an order of magnitude larger or smaller. Therefore, we write the upstream magnetic field as: 
B u = lfj,G( 1Q 3 ) B U fJi G,i9- Using this parametrization and f u = 1/2 (PN10) the various 
constraints on the maximal observed frequencies are: 

hu conf = 3 • 10 11 eV tr^-a^X^u*? 1 

hu sync = 2 • 10 8 eV tT x e£a?I% 1ia (10) 
hu IC = 3 ■ 10 5 eV r)- l Eli 4 B u ^ 19 e-V%A^Y-% 1/4 

The most robust synchrotron limit is relevant only to the GeV emission: 

'hi, \ _1 / P \ 1/4 / + \ -V4 



Implying that if the observed GeV emission is produced by synchrotron from an external 
shock in a wind environment then the acceleration mechanism must be extremely fast. This 
limit is very similar to the one obtained in the case of an ISM density profile (equation Ej), 
and is therefore general for any reasonable circum burst density profile. 

The IC limit depends on the value of the Y parameter. For the same reasons discussed 
in the ISM case, we consider here IC limits only on the X-ray emission. Unlike the ISM case 
the Y parameter of x-ray emitting electrons, Y x , depends strongly on time. In a wind density 
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profile the observed synchrotron cooing frequency, v c where most of the synchrotron energy 
is emitted, is increasing with time. As a result, Klein-Nishina effects becomes significantly 
more dominant with time, suppressing the IC co oling of X-ray emitti ng electrons at late 



time. The standard afterglow theory in a wind ( IChevalier &: Lil 120001 ) provides the value 
of v c (t) and of the Lorentz factor of X-ray emitting electrons, j x (t). Since at slow cooling 
most of the synchrotron luminosity is emitted at v c (for p < 3), KN effects are negligible for 
X-ray emitting electrons as long as 7a;Z/ c /r < m e c 2 . Thus, the time at which this inequality 
becomes an equality provides a good approximation to the time at which KN effects on Y x 
become important: 

t Yx , KN ~ 4 ■ 10 6 s Afl 7 ef A _ 2 e e ^. (12) 

This approximation assumes slow cooling and that the cooling frequency is below the X-ray 
(if the latter is not satisfied then cooling is not the limiting factor anyway, see below). It 
also ignores the effect of IC cooling on v c which can only delay the time at which KN effects 
becomes important. Thus, at t ~ 10 4 — 10 5 s, where these conditions are typically valid, 

1 /2 

Klein-Nishina effects are negligible and d _ 2 Y X is of order unity (for the same reasons 
discussed above equation [8] in the context of ISM). Thus observations of 10 keV photons 
during the first day imply: 



v 10keVy \Al, 

This result is similar to the one obtained in ISM. It implies that the conclusion that the 
observed X-ray afterglows indicate on a fast acceleration is largely independent of the circum- 
burst density profiles. The same is applicable to the conclusion that currently there is no 
strong indication for amplification of the magnetic field in the shock upstream. Note that 
this limit is valid only of the X-rays are observed to be above the cooling frequency. If X-ray 
photons are not cooling over the system dynamical time (e.g., due to a very low value of Eb) 
then the confinement limit, which require that the X-ray emitting electrons spend less time 
than the dynamical time in the upstream, is more constraining than the IC limit and should 
be used instead. 

The confinement limit is unimportant at early time, but it becomes more stringent with 
time and may become the dominant limit at very late time, t ~ 10 7 s or even earlier if 
e B <C 1CT 3 . X-ray afterglows that are observed at such late time are very rare, but they do 
exist, as we discuss in the following section. 
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4. GRB 060729 

G RB 060729 is the bu rst (at z=0.54) with the latest X-ray detection, 642 days after the 



burst (IGrupe et al. 1120101 . hereafter G10). The late time X-ray emission show a temporal 
break, from F v oc t^ 1 - 32 -"'.^) to F v oc t~^ 61 -o^\ roughly 10 6 s after the explosion. At the 
same time the X-ray spectrum varies from F u oc j^-i 1 - 18 ^ - 11 ) to F u oc z/-(°- 89 ± 011 ). This 
simultaneous temporal break and spectral hardening fits very well (within lcr) a passage 
of the cooling frequency, which increases with time, through the X-ray band (G10). This 
behavior of increasing cooling frequency is expected in a wind external medium. G10 find 
that a model of a spherical blast wave in a wind profile medium, where E = 10 54 erg (isotropic 
equivalent), A* = 0.1, Eb = 0.003 and e e = 0.1, fits the data well until t ~ 4 • 10 7 s, when a 
very sharp temporal break is observed. The spectral evolution during this late break is hard 
to constrain, due to the faintness of the signal, but it shows indications of softening. 

The origin of the late temporal break is not well determined. G10 discuss two possible 
origins - a jet break or a break in the electron distribution. They find that it is hard to 
reconcile the late break with a jet origin, although they cannot rule it out. On the other 
hand a spectral origin can provide a more consistent explanation. In that case the most 
natural source of the temporal break is inability of the shock to accelerate X-ray emitting 
particles. In that case these observations provide the first direct measurement (not only an 
upper limit) of rj. Note that according to the model of G10, the blast wave is still relativistic 
even a year after the burst, r « 5, due to the large blast wave energy and low external 
density. At late time, when the cooling frequency is above the X-ray band, the limit on 
acceleration of X-ray emitting electrons must be due to confinement. Thus if indeed the late 
break in the afterglow, at t ~ 4 • 10 7 s, is due to limited acceleration then: 

^ 100 ^- 1 (b £ o^) 1/4 ^ G ' 19 (14) 

Hence, unless the upstream field is significantly amplified by the shock precursor, acceleration 
cannot be very slow also when r « 5. Moreover, if B U fi Q ig < 0.1 then the acceleration 
must be very fast and the origin of the observed break is almost certainly due to limited 
acceleration. If B u ^c,i9 ^ 1, and the break is due to limited electron acceleration, then 
r] > 100 which is significantly larger than the value suggested by earlier X-ray (~ day) and 
GeV (~ 10 3 s) observations. This may suggest that the efficiency of particle acceleration is 
reduced when the shock approaches mildly relativistic velocities. Finally, even if the break is 
not related at all to electron acceleration then the equality in equation [T4l becomes an upper 
limit on rj. 
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5. Summary 



In this letter we examined the constraints that GRB afterglow observations poses on the 
acceleration rate, within the DSA framework, in relativistic, weakly magnetized, collisionless 
shocks in ion-electron plasma. We examine shocks that propagate into a constant density 
medium (ISM) and into a decreasing density of massive stellar wind. We consider three 
major factors that limit the acceleration in such shocks, confinement, synchrotron cooling 
and IC cooling. We find that at early times (~ 10 3 s) the best limits are set by synchrotron 
cooling of Gev emitting electrons while at intermediate times (~ 10 5 s) IC cooling of X-ray 
emitting electrons provides the best constraints. These results are independent of the circum 
burst medium density profile. At very late time (> 10 6 s) confinement may becomes the 
dominant factor in a wind environment while IC cooling remains the dominant factor that 
limits the acceleration in ISM. 

Examining available observations, the tightest limits are obtained by GeV photons that 
are observed 100 — 1000 s after the burst, if these are synchrotron photons from the external 
shock. The origin of these photons is not determined yet, but they are seen long after the 
prompt emission fades and are therefore, most likely, originate in the external shock. Various 
mode lings of the GeV emission find t hat synchrotron emission ca n explain the observations 
well jKumar k Barniol Duranlhood bold iGhisellini et al. IboiCth . If this is true, then the 
observed GeV emission require an extremely fast acceleration, at the Bohm limit or faster, 
i- e -; V ^5 !• This limit is very robust since it is almost independent of any of the shock 
parameters such as energy, density etc. On time scales of 100 — 1000 s the shock Lorentz 
factor is ~ 100 



X-ray (~ 10) keV photons are regularly observed on time scales of hours-days, where 
T ~ 10 — 50. The IC cooling of the se photons also provi des a tight limit on the acceleration 
rate: r] < 15B u ^ 5 n~ 5/8 in ISM (c.f. \u & Waxmanll2006h and r] < 30B U:llGA9 A^ /4 in a wind. 
On one hand these limits are less robust than those obtained by the GeV photons, due to the 
uncertainty in B u and the external density, but on the other hand the certainty that these 
X-ray photons are emitted by synchrotron process in the external blast wave is much higher. 
Note that if, as suggested by the GeV data (an d by the observations of other acceleration 
sites such as the Crab nebula; lAbdo et al. Il201ll ). acceleration can be as fast as rj ~ 1, then 
the available X-ray observations may be consistent with no amplific ation of the magneti c 
field in the shock upstream region (contrary to previous conclusions of lLi fc Waxmanl 120061 ). 



Finally, GRB 060729, is the burst with the longest duration X-ray afterglow observed 
to date. Its afterglow shows a sharp decline in the integrated X-ray flux 4 • 10 7 s after 
the burst (G10). This decline is most likely accompanied by a spectral softening. This 
decline can be explained naturally if the synchrotron frequency of the maximally accelerated 
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electron is crossing the X-ray band at t m 4 • 10 7 s. If this is the case than this is a direct 
measurement of 77. The afterglow light curve is consistent with a wind circum burst density 
and a cooling frequency that crosses the X-ray band at t ~ 10 6 s, implying that at later 
time the X-ray emission can be limited only by confinement. Using the fit of G10 to the 
afterglow parameters, the shock Lorentz factor at the time of the fast decline is ~ 5 and 
rj pa lOOA^ (^^Y^Bu^g, 19- If the fast decline at t > 4 • 10 7 s is not due to shock 
acceleration limit than the equality becomes an upper limit. These results suggest that the 
acceleration rate remains rather fast also at lower Lorentz factors. 

To conclude, we find that GeV and X-ray afterglow observations, provide independent 
limits on 77. The combination of these limits strongly suggest that particle acceleration is 
fast in relativistic, weakly magnetized, collisionless shocks in ion-electron plasma. Namely, 
diffusion in the shock upstream and downstream regions take place close to the Bohm limit 
at T pa 100 and it remains fast during the shock deceleration, at least up to T Rj 5. 

We thanks Rodolfo Barniol Duran, Zhuo Li and the anonymous referee for helpful 
comments. This research was partially supported by ISF grant No. 174/08 and by an ERC 
starting grant. 
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